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Abstract 


In the present work, thermo-elastic deformation analysis for the shadow mask 
assembly is performed using NASTRAN to predict beam landing shift due to the 
thermal deformation of the mask frame assembly. 

Since there exists numerous slit type apertures on the shadow mask, it is almost 
impossible to model the mask as it is. So, the shadow mask is modeled and analyzed as 
a thin plate with no apertures, introducing effective thermal conductivity and effective 
elastic constants. The effective properties are determined by subjecting a typical 
element to suitable boundary conditions and measuring appropriate quantities. The heat 
transfer in shadow mask is treated as a 2-D problem and formulation for the thermal 
analysis of a flat plate is used. The fonnulation for the plate bending analysis is used to 
calculate the thermal deformation of the mask under temperature induced strains and 
displacement boundary conditions. The three noded flat plate triangular element is 
employed for the analysis. Two test problems are solved to validate the formulation for 
2-D thermal analysis and plate bending analysis of flat plate. 

The beam landing shift of the present work (using NASTRAN) is compared with 
the results from published literature. It is observed that the temperature of shadow mask 
first increases and then stabilizes to a particular value. The beam landing shift on the 
panel first increases and then decreases to settle to a lower value. 
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Chapter 1 
Introduction 


1.1. Introduction 

A shadow mask is a thin membrane structure with numerous apertures, and is located 
behind a panel or a C.R.T. screen, so that the electron beams transiting through the 
apertures can strike the phosphors spread on the back of the screen. The thermal 
expansion of the mask frame assembly due to electron bombardment during operation 
causes a change in relative distance between the shadow mask and the phosphors. 
This results in a landing shift on the phosphor screen and, consequently deteriorates 
color purity. Predicting the temperature rise and beam landing shift is one of the most 
common requirement for designing the shadow mask type cathode ray tube (CRT). 
With increasing demands for a large size of C.R.T. with high resolution, the 
decolorization due to thermal deformation of a shadow mask has been an important 
subject of study. 

1.2. Review of literature 

Analysis of thermal deformation in a shadow mask and calculation of the 
corresponding landing shift of electron beams on the screen helps in improving the 
design of a shadow mask. A number of attempts have been made to suppress the 
thermal landing shift or mislanding on the screen. Since the development of the 
shadow mask for the color television tube [1], rimmed steel or Al-killed steel has 
been mainly used as its material. Inaba et al [1] proposed lowering the thermal 
expansion of the mask and enhancing the thermal radiation of the mask. Invar (Fe- 
36Ni) shadow mask was developed with its thermal expansion coefficient about one- 
tenth of the conventional mask. But even the Invar shadow mask CRT has a landing 
shift of the order of a few micrometer and needs to be designed precisely. ■ 
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The temperature of the mask during operation is measured by a thermocouple and 
it is assumed that the temperature of the mask is uniform over its area. Kim et al [2] 
proposed a non-contact method of measuring the temperature distribution of the 
shadow mask during operation by the radiation thermometer using an InSb photo- 
sensor to predict the correct temperature distribution in the shadow mask. 

Since the advent of shadow mask for CRT’s, bimetal springs were used to 
compensate for mask doming [3]. The bimetal system, operated by bending 
characteristics of bimetals attached to support spring, moves the shadow mask 
assembly towards the screen in such a way that the apertures follow the electron beam 
path. A comer-support system has been developed [4] that sets pins at an angle to the 
panel side wall so as to make the angle between the electron-beam path and the fiat- 
plate spring approximately equal to 90°. However, it necessitates the manufacturer to 
use a complex support mechanism to avoid disengagement of the spring and the pin 
imder mechanical loads. To compensate for the thermal mislanding [5], self-thermal 
compensator (STC), a component of the shadow mask assembly is attached to the 
frame. Still, it is difficult to control the compensation of the thermal mislanding for 
high resolution CRTs. 

Kim and Im [3] performed the finite element analysis for calculating the landing 
shift. They calculated the effective thermal conductivity and effective elastic modulus 
and developed a simplified model for estimating the surrounding temperatures. Okada 
and Ikegaki [6] simulated the landing shift by structural analysis using finite element 
method and measured data of temperature and showed that thermal deformation was 
effectively suppressed when the shape of shadow mask surface had comparatively 
large curvature on the horizontal axis and small curvature on the diagonal axis. 

One of the disadvantages of the above researches, however, is that these analysis 
procedures require measured data of shadow mask. Kim and Kim [7] analyzed the 
transient thermo-elastic deformation by FEM using the ANSYS computer program 
and calculated the beam landing shift by considering all parts inside the tube. Park et 
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al [5] simulated the landing shift by using the finite element method (FEM) 
emplo 3 dng ANSYS, thus allowing the elimination of prototyping and making the 
search for an optimal design of the shadow mask more efficient. 

In the above analysis, the shadow mask is modeled as a thin plate without slits but 
using the effective material properties [3]. A thin plate simultaneously supports 
membrane and bending actions. In classical thin plate theory, transverse shear 
deformation is neglected. However, most finite element formulation for plates and 
shells employ the first-order plate theory which accormts for the transverse shear 
deformation. When the plate is sufficiently thin and the loads are smoothly 
distributed, Reissner theory [8] is sufficiently accurate for the determination of stress 
and displacement. 

1.3. Objective of the present work 

The objective of the present work is to analyze the transient thermo-elastic 
deformation of the shadow mask assembly using the finite element method (FEM) 
through the NASTRAN computer program and calculate the beam landing shift on 
the screen. For this purpose, the shadow mask is treated as a thin plate with effective 
material properties. First, two-dimensional thermal analysis is done to determine the 
temperature distribution. After analyzing the temperature distribution, we perform the 
thermo-elastic deformation analysis using the plate bending theory. Finally, the 
landing shift of the electron beams is predicted. 

1.4. Structure of the thesis 

Chapter 2 deals with brief description of shadow mask, together with its support 
system, which compensates for the landing shifts of electron beams. The governing 
equation for heat transfer involving in-plane conduction and out of plane radiation is 
formulated and the finite element formulation is described including the nature of 
finite elements employed. Chapter 3 deals with the formulation of the transient 
thermal deformation of the mask firame assembly under temperature induced thermal 
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strains using the plate bending theory. The fourth chapter consists of the validation 
problems. The beam landing shift is predicted in the fifth chapter. The last chapter 
summarizes the work carried out and suggests the possible further scope of the 
present work. 
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Chapter 2 
Thermal Analysis 


2.1 Description of the shadow mask 

The shadow mask is used as a barrier to separate the individual color selecting 
electron beams, which are designed to strike the color phosphors [2]. Because the 
three electron beams, red, green and blue arrive at slightly different angles (from the 
three separate electron guns), see Fig. 2.1, it is possible to construct and align the 
shadow mask such that the electron beam from one gun will strike the correct 
phosphor dot, but the other two phosphors will be in shadow. This way, the intensity 
of red, green and blue can be separately controlled at each dot triad location. 


Blue Gun 



Phosphor Dots 


Red 

Phosphor 


Green 

Phosphor 


^Phosphor 


Fig. 2.1 Working of shadow mask 

The shadow mask assembly consists of a shadow mask with numerous apertures, 
frame, and support springs. The shadow mask has a great number of apertures, the 
size and shape of which vary with thickness in order to prevent dispersion of electron 
beams. The three beams are projected onto their preselected phosphors through an 
aperture in the shadow mask. The geometry of the mask is expressed by the two-step 
curvature and the aperture shape is of slit type. Because the geometry is expressed by 
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two-step curvature, Bezier surface is relied upon to express the geometry of the mask. 
The Bezier surface represents a curved surface via best approximation from vertex 
coordinates in the rectangular region, and this yields almost exact representation for 
surface with small curvatures such as shadow mask. The equation of the Bezier 
surface [9] is given as, 

m n 

(2.1a) 

<•=1 >1 

where, / (Xi,yj) are the vertices of the characteristic polyhedron that form an 
(w + 1) X (« + 1) rectangular array of points, and 
J„ i is Bernstein polynomial given as, 

=(^y(l-«)"-' ; 0<u,w<\ (2.1b) 

The shadow mask and surrounding frame are connected by spot welding at 
several points and the support springs are attached to the edge or to the comers of the 
frame. The shadow mask assembly is supported by the support springs fitted into the 
tapered stud on the C.R.T. wall. There are two types of shadow masks, when 
classified according to the way a shadow mask is connected to the frame surrounding 
the shadow mask: 

1) Bimetal system : It operates by bending characteristics of bimetals 

attached to the support spring. 

2) Comer suspension; Direct compensation system, which depends solely on 

the geometry of the support spring. 

Both of these springs are designed to compensate for the landing shift of a beam 
resulting from thermal deformation by moving the shadow mask assembly towards 
the screen in such a way that the apertures follow the electron beam path.(see Fig 2.2) 
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Panel 



Beam 


Support 

Spring 


Fig. 2 .2 C oiB^ensntion of bema luklandiug by die 
stippoit spiing 

A magnetic shield is attached to the bottom of the frame to prevent the change of 
beam trajectory due to the magnetization from the environmental magnetic field and 
to prevent the degradation of picture by stray electrons. The entire interior of CRT is 
coated with conductive coatings. The inside of the C.R.T. is in a high vacuum state 
for the electron beams to travel without being disturbed by air. Since high voltage 
difference exists between the shadow mask and electron guns, a large percentage of 
the beam energy is converted into heat energy and a little portion of beam is reflected. 
Transmitted beams through the apertures in the shadow mask do not directly 
contribute to temperature rise. The heat in the shadow mask is transferred via heat 
conduction in the in-plane directions and via heat radiation on the surface. Heat 
radiation is the more dominant heat transfer mechanism than heat conduction and the 
major parts of the heat energy ultimately flow out through radiation. The heat 
radiating into the inner walls of the C.R.T. from the surface of the shadow mask 
raises the surrounding temperature. 


2.2 Analysis procedure 

To predict beam-landing shift quantitatively, we must perform four analysis steps in 
sequence [7] : 

1. Effective material properties in the effective area of the shadow mask are 
calculated from thermal conductivity and elastic modulus in the shadow mask 
material and from the shape and layout of the apertures in the shadow mask 
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2. Transient temperature distributions of all parts inside the tube are calculated with 
the consideration of the heat radiation and the heat conduction between each part. 

3. The transient thermal deformation of the mask frame assembly is obtained from 
the temperature distribution. 

4. The beam landing shift on the phosphor screen is calculated from the relative 
distance between the panel and the shadow mask, the thermal deformation of the 
shadow mask and the position of the beam’s deflection center. 

These four steps are illustrated in Fig. 2.3 



Flow chart for Analysis 


Fig- 2.3 Flow chai-tfor analysis procedme 


2.3 Assumptions 

1. The line velocity of the electron beam with respect to the surface of the shadow 
mask is assumed to be constant regardless of position and time. 

2. The motion of electron beam is straight started from a fixed deflection center. 
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3. The heat input per unit area in the shadow mask can be obtained by taking into 
account the size of the aperture observed from the electron gun side and the 
retrace time (time during which the electron beams are not emitted from the 
electron guns) as well as the power from the electron gun. 

( 2 . 2 ) 

where, 

P : Power of electron beam, 

Amask : Total area of shadow mask, 

As : Area of aperture, 

eh : horizontal retrace time/time for scanning one line, 
ev : vertical retrace time/time for scanning the whole screen. 

PH : Horizontal pitch of the shadow mask. 

PV : Vertical pitch of the shadow mask. 

4. Since it is too complicated to consider all geometry inside the C.R.T. system, an 
simplified axisymmetric model 'as shown in Fig. 2.4 is considered [3]. We then 
replace this simplified model with an equivalent model with surrounding 
temperatures and the associated emissivities. This is also shown in Fig. 2.4. The 


^ass panel 



Fig. 2.4 4Yi!c yinitt<*tiic model for mask-frame assembly 
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equivalent emissivities and equivalent temperatures for the surrounding surface 
are determined such that the heat flux between the mask and the surroundings be 
equivalent to the heat flux in the original model. The radiation heat flux out of the 
shadow mask into the surroundings is then written as, 

=s,a(T^ -T;,) + eMT' -Ts\) (2.3) 

where, 

T ; temperature (°K) of the shadow mask, 

, £'2 : emissivities of the front and the rear surfaces respectively, 

Tsi,Ts 2 : surrounding temperatures(°K) of the front and the rear surfaces, 
respectively, 

£7 ; Stefan-Boltzmann constant (5.669 xl O'* W/m^K'*). 

5. Since there exist numerous slit type apertures on the shadow mask, it is almost 
impossible to model the shadow mask as it is. Therefore we replace the shadow 
mask by a thin plate of the same size with no apertures, introducing so called the 
effective thermal conductivity and effective elastic constants [3]. The effective 
properties are determined by subjecting a typical element (shown in Fig. 2.5) to 
suitable boundary conditions and measuring appropriate quantities. 



Fig. 1.5 Model for effective material coiustante 

a. Effective Thermal Conductivity (kx and ky) [10]; To evaluate the effective 
thermal conductivity in the x-direction, we impose a constant temperature on 
the right side and prescribe a heat flux on the left side while the top and 
bottom side are insulated. Fig. 2.6 shows these boundary conditions for 
evaluation of k^. Similar boundary conditions are applied for the evaluation of 
ky also (see Fig. 2.6). 
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Fig. 2.6 Boraidaiy conditions for evaluation of 
themal condnctmties 


Thus we have in x-direction, 

^ PV at 

Similarly, we can have in the y-direction, 

^ PH 2AT 


(2.4) 


(2.5) 


b. Effective Elastic constants (Ex,!;^, Ey»t^^) [10] : For determining the effective 

elastic constants in x direction, the lower and left sides are simply supported. 
Further a uniform displacement in x-direction is imposed on the right side 
while the constraint of the uniform displacement in the y direction is 
maintained upon the top side, (See Fig. 2.7).Young’s modulus and Poisson’s 
ratio in the x-direction are then given as. 



Fig. 2-7 Boraidaiy couditioiu; for evaluation 
of effecth'^e elastic moduli 
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II 

(2.6a) 

c 

II 

pb 

(2.6b) 

where Fx if the total reaction force in the x direction on the right side. Effective 

properties in the y direction can be found similarly as, 


F pv 

(2.7a) 

11 

(2.7b) 


6. The curvature of the shadow mask is assumed small. Therefore, it can be 
approximated as a plate. 

7. The initial temperature of the shadow mask assembly (T^ ) is taken to be 25 °C. 

2.4 Governing equation and boundary and initial conditions 

Electron beams emitted from electron gun collide with the Shadow Mask and 
generate heat energy and the heat is transferred to other parts via heat conduction and 
heat radiation. Thermal analysis of Shadow Mask is done to obtain transient 
temperature distribution of shadow mask assembly which is utilized to obtain 
deformation due to thermal forces in the shadow mask assembly and eventually 
predict beam landing shift. 

Since the thickness of the shadow mask is very small, the temperature 
gradient and heat flux in the thickness direction are almost negligible. Thus, the 
problem of heat transfer in the shadow mask can be considered as a 2-D problem 
where the incident heat flux and the radiation losses are included in the heat 
generation term. 

Consider an elementary area dA {dA = dxxdy) as shown in Fig. 2.8. Let, Q 
be the heat generation/unit area and, kx , ky be the thermal conductivities in x and y 
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direction respectively , Further, p and c be the density and specific heat 
respectively. 



Fig. 2.8 Heat exchange between faces 


Then the heat balance equation becomes, 


d ,, dT. d .. ar. ^ dT ^ 


where, r denotes the time. The, heat generation term includes the heat input into the 
bottom surface of the mask (through the electron beams) and the heat radiating out of 
the surface to the surroimding. Thus we have. 


e = k - -T‘)+s,a(T‘- T ,\ )]} (2.9) 

Since, the surrounding temperature is same for both the surfaces and the emmisivities 
of the two surfaces are also the same. 



( 2 . 10 ) 


The equation (2.9) becomes, 

Then the governing equation (2.8) becomes. 




( 2 . 11 ) 

( 2 . 12 ) 
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The boundary condition of the problem is, 

, ST ’ 

_^_ = 0 onr (2.13) 

where, F is the boundary of the shadow mask which is assumed to be insulated. The 
initial condition of the problem is, 

T = T, atr = 0 (2.14) 

2.5 Finite element formulation 

The domain is discretised into iie number of elements. Over a typical element e, the 
temperature is approximated as, 

J' = iw,(x,y)r;(r) (2.15) 

where, n„ is the number of nodes per element, TJ(r) is the time-dependent 

(unknown) nodal temperature at the local node j of the element and Nj(x,y) is the 

(known) shape function corresponding to node j. The expression (2.15) can be 
considered an approximate solution of the problem consisting of the differential 
equation (2.12) and boundary condition (2.13) if the following weighted integral is 
zero. 



(2.16) 


Here, is the area of typical element and w,- is the weighting function. Using Green’s 
lemma and the boundary condition (2.13), the first integral can be simplified. When it 
is simplified, equation (2.16) becomes. 


dw, , dT dw.- . dT 
— -k, — + — -k^ — 
dx ' dx dy dy 


3T 

lA- ^■w^q^^dA + l^■w^scr{T‘^ -T^)dA+ ^w.,pc~—dA = Q 


(2.17) 


For Galerkin’s formulation, weighting functions are given by 

(2.18) 
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Substituting equations (2.15) and (2.18) into expression (2.17), we get, 

^ dTf "" 

7=1 OT 


(2.19) 


where, 


Cy = ^pcNiNjdA , Element capacitance matrix, 

A'-’ 



dN, 
— - + 
dx 




Element conductivity matrix 

Qi = ^g,„NidA + 2 ^s&T^^^NidA. , Element heat flux vector 


( 2 . 20 ) 

( 2 . 21 ) 


( 2 . 22 ) 


Note that, while arriving at equation (2. 19) the term T* has been approximated 




( n„ 




as 

(2.23) 


J 

As a result, the matrix Ky contains 'f term where T is unknown. Thus equation 
(2. 19) becomes a nonlinear equation, which has to be solved by iteration. 


Assembling the elemental equation (2.19) over all the elements, the global 
equation becomes, 

n QT n 

7=1 M 

where n is the number of total nodes of the domain, Cy is the global capacitance 
matrix, is the global conductivity matrix, g, is the global heat flux vector and Tj 
represents the (unknown) temperature at global node j. 

Equation (2.24) is a system of ordinary differential equations. To convert it into a 
set of algebraic equations, a direct time integration technique like the forward 
difference scheme is used .For this purpose ,the time axis is divided into various time 
steps using discrete time values r*, k = 0,1,2 Let the values of nodal 
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temperatures T, at times r‘ and r“ be 7;* and T/*' respectively .Then according to 
the finite difference scheme. 


j>k+\ 

J _ ^ J 
Sr Ar 

Writing equation (2.24) at time r* 'and eliminating 
following set of algebraic equations, 


^dTj^ 


(2.25) 

using (2.25) ,we get the 


SQT/*‘=/r; ( 2 . 26 ) 

where, F/ = - LxKy)Tf + AnQ*. (2.27) 

y=i 

This set of^equation is in index form. Writing the above equation in matrix form, we 
have 


[C]{r}*^‘ = {f}* 


(2.28) 


The set of equation (2.28) is solved sequentially for the discrete times r’, 

r^,...etc. For the first time step, the right side vector is evaluated using the initial 
condition 

T] =T„ (2.29) 


2.6 Numerical integrations 

The terms of the capacitance matrix [C] and vector {F} involve siirface integrals 
which are difficult to evaluate analytically. Therefore numerical integration has to be 
applied to obtain the solution. The integrals (2.20) - (2.22) are evaluated using Gauss 
quadrature. In the present work, three-noded triangular element is used. Since the 
Gauss numerical integration formula is available only for the standard triangle, first 
the element is transformed into the standard triangle using the natural coordinates (see 
Fig 2.9). The transformation is given by. 
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(0.0) (1.0) 
stantard triangle 

Fig 2.9 Coordinate Transformation 


x = a,+a^^ + a,r]\ 

where, the constants a,- and /?, are determined from the conditions that the three 
vertices of the element get mapped onto the three vertices of the standard triangle. 
The Jacobian matrix for the transformation is given as, 

dx 

The Jacobian matrix is evaluated using the transformation (2.30) .Then the integral 
(2.20) gets transformed as 

q = I l/xN^Nj det[ ( 2 . 32 ) 

0 0 

Then, applying the Gauss numerical integration formula [1 1], we get 

c; = (2.33) 

r=I ^ 

where «c is the number of Gauss points, (^,.,7,) are the (natural) coordinates of the 
gauss points and are the corresponding weights. The integrals (2.21) and (2.22) are 
evaluated similarly: 
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•U; 1 


X-I '2 
/■=i ^ 


dN, , dNj dN. , dNj 
— -k, — - + — -k„ — - 
dx dx dy ^ dy 


+ 2s(jT^N,Nj\dGt[J]\ 




a = >*, + 2e<W,)det[y|^_^ , 

r-\ ^ 


( 2 . 34 ) 


( 2 . 35 ) 


The elements of the matrix are numerically integrated using three integration 
points located at the element midpoints of the three sides. 
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Chapter 3 

Deformation Analysis 


3.1 Introduction 

The second step in the analysis procedure is to calculate the transient thermal 
deformation of the mask frame assembly under temperature induced thermal strains 
and displacement boundary conditions. This is done by idealizing the shadow mask as 
a thin plate and using the governing equations of plate bending theory. 

3.2 Plate bending theory 
3.2.1 Basic assumptions 

Consider a plate with xy plane coinciding with the plate’s midplane and the z 
coordinate perpendicular to it and directed downwards. The fundamental 
assumptions are as follows [12]: 

1. It is assumed that the lines of material points, which are straight and normal to 
the undeformed middle surface of the element, remain straight but not 
necessarily normal to the middle surface after deformation. 

2. The strain energy corresponding to stresses perpendicular to the middle 
surface is disregarded, i.e. the stress component normal to the plate mid- 
surface is constrained to be zero in the constitutive equations. 

3 . The plate thickness remains constant during the deformation. 

Assumption (1) is equivalent to taking into account the effect of transverse shear 
deformation and it has been widely used in the context of hnear and nonlmeai 
analysis of thin plates. On the other haiid, assumptton (2) is similar to one of the plane 
stress conditions. Finally, assumption. (3) impUes that the thickness of the plate at 
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each point does not change in different deformed configurations of the plate, which 
implies zero strain along the thickness direction. 


The deformation is governed by the following three sets of governing 
equations : 

1 . Stain-displacement equations. 

2. Equilibrium equations. 

3. Constitutive equations. 


3.2.2 Strain-displacement equations 


Assumptions (1) and (3) imply that the displacement components 
arbitrary point (x, y, z) in the plate can be expressed as [13]: 

u = u^{x,y) + z(f>^{x,y) 
v = Vo(x,y) + z<zi/x,y) 


(m,v,w) of any 


(3.1) 


w = w(x,y) 


where UQ{x,y) and Vo(x,y) are the displacements of the middle surface of the plate. 
Further, and are the rotation angles of the lines normal to the undeformed 
neutral surface in the x-z and y-z plane, respectively. 


The five significant strain components may be expressed as: 


= 


Syy 


du 

dx 

dv 

dy 

1 


^xy 2 




du ^ dv 
dy dx ^ 

d\) dw 
dz dy 


(3.2a) 

(3.2b) 


(3.2c) 


(3.2d) 
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(3.2e) 


s„ = 


du dw^ 
dz dxj 


3.2.3 Equilibrium equations 

It is assumed that the body and inertial forces are negligible. Therefore, the stress 
components must satisfy the following differential equation of equilibrium, such that 
the divergence of the stress tensor must vanish. Thus, 

C7^,,.=0 (3.3) 

3.2.4 Constitutive equations 

The constitutive equation used is the generalized Hooke’s law relating the stress and 
strain components. 

where, 

: is the strain developed in the body. 

: is the initial strain developed due to differential heating. 

: Constitutive tensor. 

The expression for the constitutive tensor for isotropic material is, 

Ayw ~ ^ij^kl (3-5) 

where. A, and fi are Lame’s constants. The expression for the initial strain is, 

^oki = (3.6) 

where, 

a : Coefficient of thermal expansion. 

AT ; Temperature rise. 

5 : Kronecker’s delta. 


trncfbf sfrsfrfiirr; ?r-q-;:=r 

9Po A-- . 
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3.2.5 Boundary conditions 



Bouadaty displacements 
specified 



Bomdarytractioa 

specified 


3.1. Boundaiy conditunis 

There are two types of boundary conditions, which maybe stated as: 

1. Displacement boundary condition: Displacement is prescribed on the boxmdary. 
Thus, 


u=u 


v = v > 


on r„ 


(3.7) 



where r„ is the boundary on which displacements are specified, and {u,v,w) 

represents the specified displacement on the boundary. 

2. Traction boundary condition : Traction is prescribed on the boxmdary . Thus, 








on D 


(3.8) 


where T, is the boundary on which tractions are specified, (n ^ , ) are the 


components of unit outward normal to T, and (t„t,,t,) represent the specified 
tractions in the x , y and z directions respectively. 
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3.2.6 Variational formulation 


The principle of virtual work can be stated as, “If the displacements corresponding to 
the exact solution to the problem are perturbed by adding arbitrary virtual 
displacements, then the work done by the external forces along these virtual 
displacements equals the work done by the stresses along the corresponding virtual 
strains”. Thus, for equilibrium to .be ensured, the total potential energy must be 
stationary for variation of admissible displacements. Thus we can write:, 

SIl = {SWi-SW^) = 0 (3.9) 

where, 

n is the total potential energy for variation of admissible displacements, 

SW^ is the external work done by forces along the virtual displacements, 

SWi is the internal work done by stress along virtual strains. 


The external work done along the virtual displacements can be expressed as: 

^{txSu+tydv + tzSw)dr (3.10) 

r, 

where {du,Sv,Sw) are the components of virtual displacement. Further, the internal 
work done is given by: 

= ls(s:j - s,g jrydV (3.11) 

y 

where, V represents the domain volume, is the virtual strain and is the 
virtual initial strain. Thus, the total variation in potential energy becomes: 




- e^.j + ^)dr 




= 0 (3.12) 


3.3. Finite element formulation 
3.3.1. Matrix notation 
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Matrix Notation is used in the development of finite element equations. The 
displacement vector at any point in the domain is defined as: 

M = (3.13) 

The components of the strain tensor are represented in array form as follows: 

~ xz} O-l^) 

The components of the initial strain tensor s^ij are similarly represented as an array 


{Sg } . The components of the stress tensor (j^ are represented in array form as 
follows: 


(3.15) 

Thus, the constitutive relation (3.4) becomes, 

{o^} = [D]({s}-(s„}) (3.16) 

where the components of constitutive matrix [D] for an isotropic material may be 
written as: 


1 V 


[D] = 


E 


D 1 
0 0 

0 0 


0 

0 

1-u 

2 

0 


0 0 0 


0 

0 

0 

1-u 

2 

0 


0 

0 

0 

0 

l-v 

2 


(3.17) 


here, E is the Young’s modulus and v is poisson’s ration. 


The expression for total potential energy given by eq. (3.12) can be expressed 
using the above matrix notations, in the following form: 


sn= 


IS{{S} - {s^f[Di{s}-{ej)dV - Isiufit 


= 0 




where, 


^ f ^ t x^t y^tz 


(3.18) 

(3.19) 
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3.3.2. Flat plate element 


The flat plate element is formulated by the combination of a plane stress element and 
a plate-bending element [12], Five degrees of freedom are specified at each nodal 
point, coiresponding to three displace-ments and two rotations of the normal at the 
node. The definition of independent rotational and displacement degrees of freedom 
permits transverse shear deformation to be taken into accoimt, since rotations are not 
tied to the slope of the mid-surface.- The displacements prescribed for in-plane forces 
do not affect the bending deformations and vice versa. 


The displacement components of an arbitrary point within a flat plate element 
can be expressed as [12]: 






n 






+ z\ 


(3.20) 


where n is the number of nodes per element, Nj are the known functions of (x,y) 

known as shape functions, ,v‘ \ are the unknown displacements and 

rotations at local node i. The above equation (3.20) can be wntten as, 

{u}=[N]{uy (3.21) 


Here, the element displacement vector {«}* is given by : 

{uy = {u' ,v^ ,^y 

and, the shape function matrix is defined by: 

w.r 

{N,}\ 


where, 

{7V,}={A^, 0 0 0 Q 0 zN^ 0 . . - 

{A^2}={0 N, 0 0 zA^, 0 0 0 • •} 

{ A ^ 3}={0 0 0 0 0 0 A ^2 0 0 . . . .} 


(3.22) 


(3.23) 


(3.24a) 

(3.24b) 

(3.24c) 
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Similarly, 


(3.25) 

where, the vector 6{uY contains the nodal values of the virtual displacement vector. 

3.3.3. Finite element equations 

The domain is discretized into number of flat plate element. Then, the strain field 
is expressed in terms of nodal displacement by using eq. (3.2) and eq. (3.21) as, 

{4=[B]{»}' (3.26) 

where, [jB] is the matrix relating the strain components to the element nodal 
displacements, and is given by: 


(3-27) 


where. 


5, = 


dNi 

0 

0 

dizNi) 

0 

dx 



dx 


0 

dN; 

0 

0 

dizNi) 

dy 



dy 

dNi 

dNi 

0 

dizNi) 

d{zNi) 

dy 

dx 

dy 

dx 

0 

dNi' 

dNi 

0 

d{zNi) 

dz 

dy 


dz 

8Ni 

0 

dNi 

d{zNi) 

0 

. dz 

dx 

dz 


z=l to n (3.28) 


Similarly, 

Also, note that. 




(3.29) 

(3.30) 


Using expressions (3.25), (3.26), (3.29) and (3.30), eq.(3.18) can be written as, 

= -^{■<r[^T) = 0 P'31) 

e-1 

where, [AT]" represents the elemental stiffiiess matrix defined by, 
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( 3 . 32 ) 


[KY = \{BY[DtBW 

and, {FY represents the elemental force vector given by, 

[FT = \{Bf[D]{sMV+ |[iV]"{t'Ur (3.33) 

ye pe I J 

Here, V represents the domain of a typical element e and represents the 
elemental boundary on which the tractions are specified. 

After assembly, the above equation becomes, 

= s{uy[K]{u}-s{uy[F] = o (3.34) 

where [AT] is global stiffhess matrix and [F] is global force vector, {«} is the global 
displacement vector and S{u) is the global virtual displacement vector (or the 
variation of {«} ) . Since S{u} is arbitrary (except for the degrees of freedom lying 
on r„ ), we have 

[i:]{M} = [F] (3.35) 


The evaluation of [if]* is done using numerical integration. For this purpose, 
(x,y) coordinates are transformed to the natural coordinates (^, 77 ) using the 
transformation given by eq. (2.30). The z coordinate is transformed to ^ coordinate 
using the transformation; 


z 



(3.36) 


where t is the thickness of the plate. Thus varies varies from -1 on the bottom 
surface to +1 on the top surface. 


In the present work, three nodded flat plate element has been used for the 
deformation analysis of the mask frame assembly. The elements of the coefficient 
matrix are numerically integrated' using three integration points located at the 

midpoints of the three sides. 
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Chapter 4 
Test Problems 


4.1. Introduction 

The shadow mask assembly is analyzed using NASTRAN. To develop the familiarity 
with NASTRAN, certain test problems are solved. Results of the test problems are 
presented in this chapter. The following two test problems are solved: 

1 . 2-D heat transfer problem. 

2. Plate bending problem. 

4.2. 2-D Heat transfer problem 

An isotropic plate subjected to specified temperature on two boundaries and 
convective heat loss on the remaining boundaries [14] is analyzed using NASTRAN 
and the results are compared with those from reference [14]. This problem is chosen 
to validate formulation for the thermal analysis of a flat plate rather than for its 
practical significance. 



Fig. 4.1 Plate geometry andbomidaiy conditions 

Figure (4.1) shows the geometry and boundary conditions of the plate. Only 
quarter of the plate has been modeled due to the symmetry about the x and y axis. All 
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specified quantities refer to the quarter plate model. The geometric and material 
parameters are [14]: 

a = 400 mm. 
b = 300 mm. 
h =50E-6 W/mm^°C. 

T = 180 °C. 

To = 25 °C. 
k = 1.5E-3 W/mm°C. 

Three noded triangular element is employed for performing the thermal analysis. 
The mesh is shown in Fig. 4.2. The discretization details are: 

• Number of elements = 2400 

• Number of d.o.f = 3772. 



Fig. 4.2 Mesh generated using 3-noded triangular element 


The results obtained using NASTRAN together with the results from the 
reference [14] are presented in a tabular form in Table 4.1. Comparison is carried out 
at points 1, 2 and 3 shown in Fig. 4.1. From the Table (4.1), it is clear that the 
temperature at the three points obtained using NASTRAN are in good agreement with 


29 





the results from reference [14]. This establishes the validity of the 3 noded triangular 
element for the thermal analysis of thin plates. 


Point No. 

Temperature (°C) from 
reference [14] 

Temperature (°C) using 

NASTRAN 

1 

124.50 

123.9968 

2 1 

34.00 

39.46054 

3 

45.40 

44.30891 


Table 4.1 Comparison of temperature values 
4.3. Plate bending problem 

A plate fixed at all four edges and subjected to uniform pressure [15] is analyzed 
using NASTRAN and the results are compared using those from reference [15]. This 
problem is chosen to validate formulation for the plate bending analysis rather than 
for its practical significance. 



Pig. 4.3 Hate geometry and boundary conditions 
Figure 4.3 shows the geometry and boundary conditions of the plate. All specified 

quantities refertofull plate model. The geometric andmatetial parameters are [15]: 

L = 500 mm. 
h = 2 mm. 

E = 20pe3 N/mm^. 
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V = 0.3. 
p = 10 Pa. 


The plate is discretized using three noded triangular elements. The mesh is shown 
in Fig. 4.4. The discretization details are: 

• Number of elements = 5000 

• Number of d.o.f = 12005. 



Fig, 4.4 Mesh for plate bending analysis using 3-noded trian^e 

The results obtained from NASTRAN together with the results from reference 
[15] are presented in Table 4.2. The maximum displacement obtained from 
NASTRAN is in good agreement with the results of the reference [15]. This 
establishes the validity of the use of 3 noded triangular element for bending analysis 
of thin plates. 


Max. displacement(mm) from 
reference [15] 

Max. displacement (mm) from 

NASTRAN 

5.47E-3 

5.4E-3 


Table 4.2 Comparison of displacement values 
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Chapter 5 

Results and Discussion 


The shadow mask assembly is analyzed using NASTRAN based on the formulation 
presented in chapters 2 and 3. The three noded triangular element has been employed 
for finite element discretization. Thermo-elastic deformation analysis for the mask 
frame assembly is performed after analyzing transient temperature distribution in the 
mask. Finally, the beam landing shift is predicted. 

5.1 Problem Statement 

The shadow mask consists of a shadow mask with numerous apertures, frame and 
support springs as shown in Fig. (5.1). The shadow mask and surroimding frame are 



Fig. 5.1 Geometiy of the .shadowiiaa.sk 

connected by spot welding at several points and the support springs are attached on 
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the faces of the frame. The shadow inask assembly is supported by the support spring 
fitted into the tapered stud on the C.R.T. wall. Finite element thermo-elastic analysis 
is performed to calculate the temperature distribution and corresponding landing shift 
of electron beams on the panel due to the diermal deformation of the assembly. Since 
it is impossible to model the mask with numerous apertures, we replace the shadow 
mask by a thin plate with no apertures. Figure. (5.2) shows the shadow mask modeled 
as a thin plate. The geometric parameters, boundary conditions and material 
properties are ; 

1. Geometric parameters [3]: 

Length of mask, A = 404.8 mm. 

Length of portion with apertures, a =387.92 mm. 

Width of mask , B = 306.0 mm. 

Width of portion with apertures, b = 298.12 mm. 


Thickness of mask , h = 0.2 mm. 



Fifi. 5.2 Thin plate model for shadow mask 


2. Boundary Conditions 

. Input heat flux qin [3]= 210. 1 W/mml 
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• Thermal boundary conditions: 

Insulated at edges. 

• Geometric boundary conditions 

Clamped at the edges, u = 0 , v = 0, w=0. 

3. Material Properties : 

The effective material properties for the mask without apertures and the mask with 
apertures are shown in Table (5.1) [3] . 


Properties 

Mask 

( without apertures) 

Mask 

(with apertures ) 

Young’s modulus 
(N/mm^) 

2.2E5 

1.6E5 

Poisson’s ratio 

0.25 

0.26 

Specific heat 
(J/kg°C) 

475 

475 

Conductivity 
(W/mm °C) 

5.45E-2 

3.055E-2 

Co-eff of thermal 
expansion ( / °C) 

1.15E-5 

1.15E-5 

Emmisivities 

0.55 

0.87 


Table 5.1 Material properties 


The geometry of the plate is discretized using the three noded triangular element 
for both the transient thermal analysis as well as deformation analysis. The 
discretization details are: 

Number of elements = 3584. 

Number of d.o.f = 5613. 

5.2 Transient thermal analysis 

Transient thermal analysis is performed on the plate model to determine the transient 
temperature distribution in the shadow mask. Figures (5.3 - 5.7) show the 
temperature distribution on the shadow mask after 50 seconds, 130 seconds, 250 
seconds, 410 seconds and 490 seconds of operation respectively. A steady state is 
reached approximately at 490 seconds. The maximum temperature occurs at the 
center of the mask. Figure 5.8 shows the variation of maximum temperature with time 
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5.8S4^: 



Fig. 5.3 Temperatiu’e <listiibution (at t = 50 sec) 



Fig. 5.4 Temperatiue clistiibxitioii {att= 130 sec) 




y 

z" 


Fig. ff.S Temperature distiibution (at t = 250 sec) 



Fig. 5.6 Temperatiue distiibutioii (at t — 410 sec) 






Fig. 5.8 Temperatiu'e at center of mask using NASTEAN 
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Fig 5.9 Compaiisoii of teinijeratiue (U.stiibiitiou'vritii reference {3] 

Figure 5.9 compares the results from the present work with the result of reference 
[3]. From Fig. 5.9 it can be inferred that the trend of temperature distribution for the 
present work is in agreement with reference [3]. 

5.2 Thermal deformation analysis 

Thermal deformation analysis is performed on the plate model for the shadow mask 
to estimate the deformation of the mask corresponding to the thermal strain induced 
in the mask. Figures (5.10 — 5.14) show the contours of beam landing shift on shadow 
mask after 50 seconds, 130 seconds, 250 seconds, 410 seconds and 490 seconds of 
operation respectively. Beam landing shift is the magnitude of resultant displacement. 
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Fig. 5.15 Beam lamlhig sliift at position.^ 1 (18$, 133) and 2( 123, 0) 
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Figure 5.15 shows the variation of beam landing shift with load cases. The load cases 
correspond to the time steps taken for the transient thermal analysis of the mask. The 
peak value of beam landing shift (3.48 E-3 mm) occurs at 90 seconds, and thereafter 
the beam landing shift settles to a lower value as the temperature difference between 
the center of mask and comer of mask begins to decrease. The present result differs 
appreciably from that of reference [3]. This discrepancy in the result is due to the 
difference in the modeling. The shadow mask here has been modeled as a flat plate 
rather than as a curved surface along with the frame and support springs which will 
present the actual behavior of the shadow mask response to the thermal deformation. 
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Chapter 6 

Conclusions and Scope for Future Work 


6.1 Conclusions 

The shadow mask assembly is modeled as a thin plate for predicting beam landing 
shift due to thermal deformation of the mask. The heat transfer in shadow mask is 
treated as a 2-D problem and formulation for the thermal analysis of a flat plate is 
used. The formulation for the plate bending analysis is used to calculate the thermal 
deformation of the mask under temperature induced strains and displacement 
boundary conditions. Two test problems are solved to validate the formulation for 2- 
D thermal analysis and plate bending analysis of flat plate. 

The three noded triangular element has been employed for finite element 
discretization to perform thermo-elastic analysis using NASTRAN. The trend for the 
temperature distribution in the shadow mask using NASTRAN is in good agreement 
with the temperature distribution from literature. The trend for the beam landing shift 
is also in agreement with landing shift fi'om the literature, thou^ the values differ. 
This discrepancy in the result is due to the difference in the modeling of the shadow 
mask. This shows that the assumption of modeling the shadow mask as a thin plate 
can be used to predict the trends of temperature distribution on tiie mask and the 
beam landing shift but not their values. 

6.2 Scope for Future Work 

The present work can be extended in the following directions; 

1. The shadow mask can be modeled as a shell and the governing equations for 
heat transfer and deformation for shell can be used. 
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2. The whole shadow mask assembly comprising of shadow mask, frame and 
clips should be analyzed together to precisely predict the landing shift. 
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